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We show the result of the lattice gauge field tensor which is derived from the classical continuum 
limit of the overlap Dirac operator. By analogous construction, it was recently proposed that the 
gauge action can be obtained from the overlap operator as well which is proportional to TraD ov . 
We discuss how to carry out Monte Carlo simulations of such a gauge action together with the 
overlap fermion. 
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1. Introduction 

Lattice gauge actions are usually constructed from the link variable U. Similarly, the gauge 
field tensor is constructed from the clover-leaf like plaquettes. 

The advent of Neuberger's overlap operator [jl]] has the implication in a new direction. The 
overlap Dirac operator is development to solve the fermion chirality problem in lattice QCD [Q]. It 
does not have 0(a) error for the fermion operators. The small 0(m 2 a 2 ) errors and non-perturbative 
renormalization via current algebra relations and Ward identities make it a desirable fermion for- 
mulation for both light and heavy quarks [Bj. An index theorem was formulated by Hasenfratz, 
Laliena, and Niedermayer [|j] based on the operator at finite cut-off. It is then pointed out by 
Liischer that the anomalous behavior of the fermion partition function under a flavor-singlet trans- 
formation is expressed by the index of the the overlap Dirac operator arising from the Jacobian, 
providing a clear understanding of the exact index theorem [||, in the path-integral formalism. 
Following developments have seen the explicit derivation of the local lattice topological charge in 
terms of the local overlap operator via weak coupling expansion by Kikukawa and Yamada [g], ex- 
plicit calculations without gauge coupling expansion by Adams [f7|], Fujikawa [j|] and Suzuki [g], 
i.e. 

q L (x) = tr c ,y 5 (l - l/2aD m (x,x)) (1.1) 

where D ov is the overlap operator and the trace is over the color and Dirac spin indices. The lattice 
topological charge operator #z,(jc) so defined approaches the topological charge density q(x) at the 
continuum limit, 

q L (x)-^a 4 q(x) + 0(a 6 ). (1.2) 
This formulation of the topological charge operator has not only been used to check the Atiya- 



Singer index theorem at finite lattice cut-off [10], the topological susceptibility [11], but has also 



been adopted to study the local topological structure of the vacuum [|l2|, 13, 14, 15, T(| 17, |T|]. 
It is with this operator that the sub-dimensional coherent sign structures have been discovered in 
4-D QCD 00,110] and 2-D CP(N-l) model |g9||. It is argued that the chiral filtering of 



the overlap fermion is responsible for optimally filtering out the ultraviolet fluctuation to allow the 



structures to be revealed [20]. Indeed, other conventional topological charge operator constructed 
from the gauge links were used, but could not decipher the curvilinear structure observed with 
the overlap operator |Jl9[]. This leads to the possibility that gauge field tensor and gauge action 
derived from the overlap operator might be good alternatives to those from the gauge links directly. 
Recently, I. Horvath [ pl| ] proposed a formulation of lattice QCD wherein all elements of the theory 
(gauge action, fermionic action, theta-term, and other operators) are constructed from a single 
object, namely the lattice Dirac operator D ov . In this talk, I will present the results of the derivation 
of the gauge field tensor as the classical continuum limit from the overlap operator. The detailed 



derivation [22] and the numerical evaluation [23] are under preparation and will be posted on 



the arXiv soon. I will also discuss how to simulate such an action with the Hybrid Monte Carlo 
algorithm. 

2. Gauge Field Tensor 

To begin with, we note that a gauge covariant operator which is a functional of U satisfies the 
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condition 



0(g- l Ug)=g- l O(U)g, 



(2.1) 



where g is the local gauge transformation. It is obvious that the local operator with the color trace 
tr c O(U)(x,x) is gauge invariant 



tr c O{g- l Ug)(x,x) =tv c g- l {x)0{U)(x,x)g(x) =tr c O(U)(x,x). 



(2.2) 



Since the overlap operator, being a Dirac fermion operator, is gauge covariant and is not ultra- 
local, it is expected that the classical continuum limit of the trace in both the color and spin indices 
tr cj rD 0V (jc,x) will be the lowest dimensional local gauge operator which reflects the Lorentz struc- 
ture of the gamma matrix T. Thus, it is not surprising that tr cv YsD ov (x,x) gives the local topological 
charge density at the continuum limit. Therefore, one expects that [ gjj ] 



tr cs (D ov (x,x) - D° ov (x,x)) a 4 tr c F 2 v + 0(a 6 ), 



(2.3) 



where D° , is the non-interacting overlap operator with gauge link U set to unity. This has been 



verified with constant field and with numerical evaluation [ 23 ] 



In addition to gauge invariant operators, one can obtain gauge covariant operators as well. 
Since D ov (x,x) is gauge covariant and a Lorentz scalar, one expects that tr s a^ v D ov (x,x) with the 
spin index traced over will result in a lowest dimensional gauge covariant operator with the /iv 



indices in the classical continuum limit which is just the gauge field tensor [21]. In other words, 



tr s a^ v D 0V (x,x) ^ a 2 F^ v + 0(a 4 ). 



(2.4) 



We should note that the possibility of obtaining the lattice gauge field tensor from the overlap 



operator by expanding in a was first pointed out by Niedermayer Q24Q. Also, it was suggested 



by Gattringer [ [250 that the field tensor can be defined through the square of various lattice Dirac 
operators, i.e. tr s Y^j v D 2 (x,y) with a weighted sum over y. 

We have derived the classical limit following Adams [j7|] and Suzuki's method [^]. While 



details of the derivation will be given elsewhere [22], we shall present the results here. 



tTgffuvODovC*,*) — 

d 4 k -2 



n (27T) 4 (ZtZ)3/2 



{m + rYfcx - \)c^c v + 2rc^s 



where = sink^,c^ 



cosk^ and 



in 



+ rE(c M -l) 



a 2 F^ v (x) + 0(a 4 



(2.5) 



(2.6) 



For the case where r = l,m = 1.368 (which corresponds to Wilson fc = 0.19 in the kernel), 
the proportional constant is —0.08156 for the overlap fermion [E5]. This has been numerically ver- 



ified [23]. One can try to use this operator to calculate glue properties, such as the glue momentum 
and angular momentum fractions in the nucleon to see if they are better than the glue operators 
constructed from the link variables as is in the case of the local topological charge. 
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3. Monte Carlo Simulation 



Take the proportional constant for the classical continuum limit in Eq. ( |2.3[ ) to be c which has 
been evaluated [ [23| ] as a function of the negative mass parameter p in D w , the kernel of the overlap 
operator D ov . Then the lattice gauge action can be written as 

! Tr(D ov -D° ov ), (3.1) 



leg 



where Tr denotes the trace over all indices, e.g. color, spin, and space-time. For a given lattice, 
the trace over the free quark overlap operator, i.e. TrD m is a constant which has no effect on the 
Markov process to obtain equilibrium gauge configurations, we shall drop it. Noting that D ov obeys 
75 hermiticity and satisfies Ginsparg-Wilson relation, we have 

Sg = 2^ TrDov = {Dm + D - } = 4^2 TrD - D - (3 " 2) 

Now the lattice QCD partition function is 

Z = J ®U d%dy f e- S * +1: % V^K))^ (3 .3) 

where D ov (mf) is the overlap operator for a quark with mass mj 

D ov {m f ) = pD ov + m f {\ - -D m ), (3.4) 

where p is the negative mass parameter in D\y- Since e TrM = e Trlne " = dete M , the gauge part of 



the partition function in Eq. ( |3.3| ) can be written in terms of a fictitious fermion field \jf g so that 



p _ /V j- 

Z = mdy g d ¥g dxjr f d ¥f e^ e Acs ~ )%+E/iiW(o»vK))v/. (3 . 5) 



After integration of the fermion fields, it can be written as 

N f 

n 



r 1 t Nf 

J ®Ufet{e~^ D ° vDov )]\&et(D ov (m f )). (3.6) 



We see that the gauge action plays the role of a UV-filtering for the fermion action (note that 
c > for the range of parameters for the Wilson kernel in the overlap operator with r = 1 and 



2 > m > 1 [|23[]). The efficiency of an UV-filtered fermion determinant has been studied by Duncan, 
Eichten, and Thacker [^] and by Borici [|27|]. 

One way to carry out Monte Carlo simulation is to use pseudofermions to simulate the deter- 
minant. For example, one can equally split the gauge determinant and attach them to the fermion 
determinants of different flavor. In terms of the pseudofermions, it is 



— — jD' ov D G v 



9Udtfdtye-' L &*r e *' D -K)^/. (3 . 7) 
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We shall discuss two ways to approximate the pseudofermion action. Since D oc is normal, i.e. 
[Dl v ,D 0V ] = 0, one can write 

f f e^^D-J (m f » f = r f ~ l " K n,J nn l H— l, "-O l , (3.8) 

consider the Chebyshev polynomial approximation to degree M 

J_nt n M 

I 

i'=l 



The range of eigenvalues of D^ V D 0V is from to 4p 2 . If = ' 2 is about unity or less, one can 

ACLVfg 



°" " -t^CiiDlpJ. (3.9) 



Alternatively, one can perform a Chebyshev rational polynomial approximation for the operator 

e™f* tDmDm D^{mf) to degree N for the Rational Hybrid Monte Carlo algorithm (RHMC) [§f]. 
For the case of 2 + 1 flavors, the pseudofermion action for the 2 degenerate flavors can be approxi- 
mated by 

re^ DlD -(DlD 0V (m f ))- l <t> ~ 0* f , - 0, (3.10) 

i=l D ] ov D ov + bi 

and the single flavor one approximated by 

r e ^P Dl ^{DlD m {m f ))- 1 ^ ~ 0* £ Cf f (3.11) 

i=l D ov D ov + di 

In this case, the forces in the equation of motion in HMC come from the effective pseud- 
ofermion actions in Eqs. ( 3.10) ) and ( 3.11 ) which represent the combined gauge and fermion forces. 



If one uses a multi-mass algorithm for inversion and the coefficients bi and d, are not smaller 
than rrrj, the overhead of incorporating the gauge action in RHMC is negligible compared to the 
ordinary 2+1 flavor simulation in HMC with the same inverter. 

One can of course accelerate the Rational Hybrid Monte Carlo algorithm (RHMC) by splitting 



the determinant into fractional flavors [g^, [30p with more pseudofermion fields and improve the 
overall efficiency as shown by Clark and Kennedy [f!^]. 

In summary, we have discussed possible Monte Carlo simulations of the lattice gauge action 
from the trace of the overlap operator, i.e. TrD ov together with the overlap fermion action in the 
context of HMC. By virtue of the fact that the overlap operator is exponentially local, the gauge 
action so defined is expected to behave like a chirally smeared action. Furthermore, the integrand 
of the gauge part of the partition function, written in terms of a determinant, appears to be an UV- 
filter for the fermion determinant. We should note that, similar to the overlap fermion action, this 
gauge action is not reflection positive. Also presented is our derived result of the lattice gauge field 
tensor as the classical continuum limit of tr s a^ v D ov (x,x). This can be used to calculate glue matrix 
elements in the hadrons and possibly glueballs. 

This work is partially supported by DOE grant DE-FG05-84ER40154. We wish to thank 
I. Horvath, A. Alexandur and A. Kennedy for stimulating discussions. The author also wish to 
thank the hospitality of the Institute of Physics, Academia Sinica, Taipei, Taiwan, National e- 
Science Center, Edinburgh and DESY-Zeuthen, Germany where the talk is prepared and the pro- 
ceedings is written up during the visit of these places. 
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